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Abstract. We study the spreading of a disease on top of structured scale- free networks recently introduced. 
By means of numerical simulations we analyze the SIS and the SIR models. Our results show that when 
the connectivity fluctuations of the network are unbounded whether the epidemic threshold exists strongly 
depends on the initial density of infected individuals and the type of epidemiological model considered. 
Analytical arguments are provided in order to account for the observed behavior. We conclude that the 
peculiar topological features of this network and the absence of small-world properties determine the 
dynamics of epidemic spreading. 
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1 Introduction 

During the last years, there has been a burst of activ- 
ity in the study of complex networks It has been 
shown that many social and natural systems |0,0,^,^,Q,pl 
^,[uj[ll) can be represented as a graph where nodes repre- 
sent individuals or agents and links stand for the physical 
interactions among them. Surprisingly, many of these net- 
works share some important topological features such as 
' small-world (SW) properties jl2j and scale-free (SF) de- 
gree distributions p3| , where the degree or connectivity 
k of a node is the number of neighbors to which it is 
linked. Networks displaying scale- free degree distributions 
are very interesting not only for their relative abundance 
in Nature but also because of their peculiar statistical 
properties. In particular, the unbounded fluctuations of 

■ the connectivity distribution P(k) ~ fc~ 7 (7 < 3) seem to 
be a blueprint of all real-world networks studied so far jj], 

!• 

The SW and SF properties have considerable impact 

■ on the processes running on top of complex networks. The 
effects of the complex topological features of SF networks 
on the dynamics of epidemic disease spreading is per- 
haps one of the most interesting outcomes in the study 
of complex networks. Specifically, the understanding of 
spreading phenomena in these networks can shed light on 
a large number of practical problems, ranging from com- 
puter virus infections (l4|,|l5| to epidemiology [^6i|l7|l . For 
instance, it has been recognized only recently [|18| , |1S[| that 
random SF networks are completely prone to epidemic 
spreading allowing the onset of large epidemics whatever 



the spreading rate of the infection. This radical change 
with respect to the way a disease is spread in a regular 
structure is rooted in the diverging connectivity fluctu- 
ations of SF networks with 7 < 3. Similarly, random SF 
networks have been s hown to e xhibit extremely robustness 
to random damages 



21 



e-mail: yamir@ictp.trieste.it 



On the other hand, real networks are also character- 
ized by degree correlations that might play a fundamental 
role in the functional properties of networked systems and 
on processes running on correlated networks. The study 
of networks with degree correlations and the extension of 
previous results obtained for random networks to corre- 
lated networks is quite recent g|^,|^|^,|2^,|8]]. One of 
the models aimed to account for degree correlations was 
proposed recently by Klemm and Egufluz p9| (henceforth 
referred to as deactivation model) where nodes can be de- 
activated with probability inversely proportional to their 
connectivity. The model turns out to be more rich than 
initially thought with a power law connectivity distribu- 
tion P{k) ~ fc~ 7 but with 7 laying in a quite wide inter- 
val 2 < 7 < 4 p0[ . Additionally, it has been shown that 
the topology of the network is essentially linear being the 
graph a collection of stars of diverse degree connected as a 
chain by a number of local links, i.e., the grap h lacks SW 
properties. While a threshold may exist |3l| due to the 
peculiar topological properties of the network that makes 
the epidemic spreading to be dominated by the diffusion 
of the disease on a linear chain J3fj, the star-like graphs 
connected as a chain might also lead to new effects in the 
dynamics of epidemic spreading. 

In this paper, we study in detail by large scale numer- 
ical simulations two paradigmatic epidemiological models, 
namely, the Susceptible-Infected-Susceptible (SIS) and the 
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Susceptible-Infected- Removed (SIR) models on top of net- 
works generated using the deactivation model. We found 
that the existence or not of an epidemic threshold for 
7 < 3 depends on the initial density of infected individ- 
uals while for values of 7 > 3 the epidemic threshold is 
recovered and it is universal, i.e., it does not depend on 
the initial density of infected individuals. We also provide 
analytical arguments in order to explain our numerical 
findings. Moreover, we show that the linear character of 
the graph determines the existence of a trivial threshold 
for the SIR model that can be mapped to a bond perco- 
lation problem in one dimension. 



2 Deactivation model 

Let us first summarize the main properties of the deacti- 
vation model introduced by Klemm and Egufluz |p9|] . The 
growing dynamics of the network is defined in the follow- 
ing way: One starts from a fully connected graph of m 
nodes that are set active. At successive time steps nodes 
are added one by one following the recipe: (i) A new node 
is connected to all active nodes in the network; (it) one 
of the active nodes is selected for deactivation with prob- 
ability 



(1) 



where the sum in Eq. (|l|) runs over the set of active nodes 
A, a is a model parameter, and kf 1 denotes the in-degree 
of the z-th node. Finally, (Hi) the new node just added is 
set active. 

Recently, it has been shown |3(| that most of the topo- 
logical properties (in particular, the degree distribution) 
of the network generated using the above rules are very 
sensitive to the order in which steps (ii) and (Hi) are per- 
formed. Thus we shall discriminate in the following two 
cases 

— Model A: {ii) is performed before {Hi) . 

— Model B: (Hi) is performed after {ii). 

The deactivation model is usually run taking a = m. 
This makes the deactivation probability inversely propor- 
tional to the total connectivity of the nodes k = m + k ln . 
The connectivity distribution can be analytically obtained 
for the limiting cases of lowest |Q and largest m [|T] for 
each model resulting in a power law that reads as, 



P{k) ~ fc~ 7 , 



(2) 



where the exact value of 7 depends on the model consid- 
ered and the value of m such that. 



model A with a = m 
model B with a = m 



3 < 7 < 4 
2 < 7 < 3. 



Noticeable and relevant for the physical processes run 
on top of this model, the dynamics of the deactivation 
model allows exponents 7 of the connectivity distribution 
that can lead to an unbounded second moment (Ik 2 ) — > 00 



for model B) or to finite connectivity fluctuations ((k 2 ) is 
finite for model A). More important, however, is the fact 
that the generated networks lack SW properties. Small- 
world properties refer to the fact that in many complex 
networks one can go from one node to any other node of 
the system visiting a very small number of intermediate 
nodes. The minimum number of such intermediate nodes 
for each pair (i,j) of nodes in the network is called the 
minimum path length between i and j, and the diameter 
of the network is defined as the largest among the shortest 
paths between any two nodes in the network. For networks 
displaying SW properties, the graph diameter grows, at 
most, logarithmically with the system size N. 

The networks generated with the deactivation model 
consist of a collection of star-like graphs connected as a 
chain, resembling the topology of a one dimensional lat- 
tice. In fact, irrespective of the version considered, both 
the diameter and the average minimum path distance scale 
linearly with the network size N |so| . In particular, the 
mean square displacement of a random walker on the de- 
activation model scales with time almost as for a one di- 
mensional lattice qualifying for a slightly subdiffusive be- 
havior |30|. Hence, the existence of an epidemic threshold 
in the deactivation model is not surprising. However, at 
the same time, the peculiar characteristics of the graph 
and, in particular, the presence of stars with very large 
connectivities introduces new and interesting effects in the 
dynamics of disease spreading. As we shall discuss in the 
following section, whether or not a finite prevalence sets 
in for the SIS model would depend on the initial density 
of infected individuals. 



3 The SIS Model 

In the standard SIS epidemiological model jl6| , each node 
of the network represents an individual and each link is 
the connection along which the individuals interact and 
the epidemic can be transmitted. In this model, the indi- 
viduals can exist in two possible states, namely, suscep- 
tible or healthy and infected. Individuals removal due to 
death or acquired immunization is not allowed and thus 
they stochastically move through the cycle susceptible — > 
infected — > susceptible. The disease transmission is de- 
scribed in an effective way. At each time step, each sus- 
ceptible node connected to one or more infected nodes 
gets the infection with probability A, while the infected 
nodes are cured and become again susceptible at a rate 
S (henceforth set to 1 without lack of generality). In the 
SIS model whenever the epidemic pervades the system, it 
gets into an endemic state with a stationary value for the 
density of infected nodes that acts as the order parameter 
of the model. 

Results on regular structures point out that there is a 
threshold below which the system does not reach such a 
stationary state and the infection dies out p6| , |TS| ] , More- 
over, the value of the prevalence at the stationary state 
and its own existence does not depend on the initial den- 
sity of infected individuals. One can start with an initial 
density po — 1/N of infected individuals or from, let's say, 



Yamir Moreno, Alexei Vazquez: Disease Spreading in Structured Scale-Free Networks 



3 



half of the lattice infected (in general, po = const): The 
prevalence self-organizes into a stationary state in both 
cases. It is worth noting that both situations are feasible 
in practice and thus physically relevant. 

On the other hand, the behavior of the SIS model on 
random SF networks is radically different |l8],[l]| when 
7 < 3. In this case, there is no epidemic threshold in 
the thermodynamic limit and the networks are completely 
prone to the spreading of the disease. This drastic change 
of behavior is due to the diverging connectivity fluctua- 
tions of SF networks. Moreover, it has also been shown 
that for finite-size system an effective threshold is recov- 
ered, but its value is significantly overestimated |3^] . 

In order to study the SIS model in deactivation model 
networks, we first generate networks using the algorithm 
described in the preceding section for both settings of the 
model A and B with a = m = 3. Then we let the system 
evolve according to the SIS dynamics. Large scale numer- 
ical simulations were performed in networks of size up to 
N = 10 6 , averaging over at least 100 different realizations 
on at least 10 different realizations of the network. Ini- 
tially, a fraction of nodes po was infected and we let the 
system relax into the steady state where the prevalence p 
attains its stationary value. Two initial conditions for the 
density of infected individuals were considered. In the first 
case, we start from a single infected individual, po = 1/AT, 
and set the observation time t m to be 10 6 time steps, that 
is, the state of each node is updated up to t m times if 
there is at least one infected individual. Then, we repeat 
the process but starting with homogeneous initial condi- 
tions. This is achieved by infecting a finite fraction of the 
network po = const, and monitoring as before whether a 
stationary state sets in with its corresponding prevalence 
for the same observation time t m . It is worth recalling that 
the prevalence in the stationary state is computed as the 
average over all surviving trials in both cases. 

The results obtained are depicted in Fig. |l| where the 
steady state prevalence has been drawn as a function of the 
spreading rate A. As expected for model A the epidemic 
threshold would appear again for both initial conditions 
since in this case together with the lack of SW proper- 
ties and the linearity of the network, the second moment 
(fc 2 ) is finite with 3 < 7 < 4. In this case, the phase 
diagram is the same regardless of the initial density of in- 
fected individuals po- It is worth noting that whenever an 
epidemic threshold exists, its value seems to be the same 
(within statistical errors) for both versions of the deacti- 
vation model, depending only on the average connectivity 
(k). This implies in its turn that the dynamics of the infec- 
tion spreading is mainly determined by the linear star-like 
structure of the graphs. 

Two radically different behaviors can be noted for the 
deactivation model B depending on the initial condition. 
In this case the exponent of the power law decay of the 
degree distribution is in the interval 2 < 7 < 3. In par- 
ticular a best fit analysis for networks of size N = 10 7 
yields an exponent 7 = 2.65 ± 0.05 for m = 3. When 
the system evolves from a single infected individual, an 
epidemic threshold A c appears, below which the epidemic 
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Fig. 1. Phase diagram for the deactivation model A (top) 
and B (bottom). The density of infected nodes p for two initial 
conditions is plotted as a function of the spreading rate A. The 
size of the system is N = 10 6 and the average connectivity was 
set to (k) — 6 (m = 3). For model A the epidemic threshold 
is independent of po and is equal to A c = 0.16(1). For model 
B the epidemic threshold depends on po. For an initial state 
with a single infected individual A c = 0.14(1), whereas for ho- 
mogeneous initial conditions A c — » in the thermodynamic 
limit. 



cannot pervade the network. In this case the fact that in 
the thermodynamic limit the connectivity fluctuations are 
diverging is canceled by the effects introduced by the lin- 
ear topological nature of the network and the lack of SW 
properties. As advanced in Ref. |3(|, the SIS dynamics 
can be reduced in a coarse grained picture to the diffu- 
sion of the disease on a linear chain. Thus, the classical 
picture for regular structures with an epidemic threshold 
is recovered. 

For homogeneous initial conditions, the phase diagram 
is completely different. In this case, the stars with high 
connectivity can let the epidemic survive and a final sta- 
tionary state is reached by a long power-law decay. In 
order to rule out the presence of finite size effects hiding 
an abrupt transition we have checked the behavior of the 
stationary prevalence for several system sizes. Note that 
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Fig. 2. Prevalence p for the SIS model in scale-free correlated 
networks generated using deactivation model B. The average 
connectivity of the network is in this case (fc) — 6 (m = 3). 
The straight line is a fit to the form p(A) ~ A* 9 , with j3 ~ 4.5. 



larger system sizes are needed in order to correctly de- 
pict the prevalence behavior for A <C 1, since finite size 
networks induce finite size corrections to the zero thresh- 
old fl|,(L9|. The results obtained are plotted in Fig. |. 
The straight line is a fit to the form p(X) ~ A* 3 , with 
[3 ~ 4.5, showing that the infection prevalence is assum- 
ing a finite stationary value for all values A > 0. This nu- 
merical evidence confirms the previous picture obtained 
for uncorrelated scale- free networks, namely, the absence 
of an epidemic threshold in the thermodynamic limit. It 
is interesting to recall that using a dynamical mean field 
approach for uncorrelated networks with scale-free con- 
nectivity distributions given by P{k) ~ fc~ 7 , the same 
functional form p{\) ~ A' 3 is predicted when 2 < 7 < 3, 
although in this case satisfying the relation (3 = 1/(3 — 7) 
Hence, it seems that correlations, while preserving 
the general functional form, make the prevalence to decay 
more faster than it does in the uncorrelated case. This 
in turn may explain why in finite systems the epidemic 
threshold is smaller than its counterpart in uncorrelated 
networks, a fact also observed in random network models 
Ipij . Below, we provide an analytical argument that help 
understand the origin of these differences by analyzing the 
SIS model on an ensemble of disconnected stars. 



3.1 SIS model on an ensemble of disconnected stars 

Let us study a graph made up of an ensemble of discon- 
nected stars. To be more precise, we consider a graph 
where all nodes with degree k > 1 are connected to nodes 
with degree k = 1. If P(k) (k > 1) is the degree distribu- 
tion then P(l) = J2k>i kP{k). Since the stars are isolated 
the stationary prevalence is given by 



P : 



(3) 



where n(k) is the stationary number of infected nodes on 
a star of degree k and Vk is the probability that at t = 
there was at least one infected node in a star with k leaves. 

Now, to compute n(k) we focus on the SIS dynamics on 
a single star. The spreading takes place from the central 
node to the leaves, then from the leaves to the central node 
and so on. For stars with a large number of leaves (k 3> 1) 
we can approximate the average number of infected nodes 
by the number of infected leaves. Hence, the number of 
infected nodes at step t + 1 is 



n(k,t+l)= l-(l-A) 



Afc 



(4) 



The first factor in the right hand side is the probability 
that the central node gets infected, receiving the disease 
from at least one leaf. The second factor is the number 
of infected leaves given that the center was infected. This 
linear map has always the trivial solution n(fc,i) = 0. To 
investigate its stability we assume n(fc,i) <C 1 resulting 



n(k,t+ 1) = Afc In- 



1 



1 - A 

Thus, there is a critical star size 



n(fc, t) 



fe c (A) = A In 



1 



1 - A 



(5) 



(6) 



such that for fc < fc c the prefactor in the right hand side 
of Eq. (||) is smaller than 1 and, therefore, n{k) = 0. On 
the contrary, for k > fc c there is an exponential growth 
of n(k 7 t) indicating that the solution n(k,t) = is not 
stable, i.e. n{k) > 0. Moreover, the critical degree depends 
on A, with the limiting cases fc c (0) = 00 and fc c (l) = 0. 
Going back to Eq. (||) we obtain that 



P(k)n(k)v(k) 



(7) 



Let us now distinguish between two different initial con- 
ditions considered before. In one case we start from only 
one infected node. The probability that one node of a star 
of k leaves is infected is 



;W(fc) = 



k + 1 

N 



The substitution of this expression in Eq. (Q) yields 



k>k c 



N 



(8) 



(9) 



Now, for 2 < 7 < 3 the second moment grows at most 
as ^fc 2 ) = 0(7V 3 ~ 7 ). Hence, in the thermodynamic limit 
the prevalence is asymptotically zero. However, a different 
result is obtained if we infect a finite fraction po of the 
nodes. In this case the probability that at least one node 
of a star with fc leaves is infected is given by 



fc>i 



,(2) 



(fc) = 1 - (1 - P o) k+1 > po 



(10) 
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From this expression and Eq. (Q) it follows that 

p (2) ^p°Y. p ( k M k ) ■ (11) 

Thus, there is a finite prevalence in the stationary state. 

These two opposite scenarios are the same observed in 
the numerical simulations of the preceding section. When 
there is a finite fraction of infected nodes at t = 0, a finite 
prevalence is obtained. On the contrary, when the initial 
infection is concentrated in only one node, the station- 
ary state reflects the topological nature of the graph. For 
the ensemble of disconnected stars we obtain a station- 
ary state with no infected node, corresponding to the SIS 
model in a zero-dimensional system. On the other hand, 
for the deactivation model, we get a phase transition at 
a finite infection rate as expected for a one-dimensional 
system [p3[. 



4 SIR model 

The SIR model assumes that individuals can exist in three 
classes: susceptible, infected and removed. The main dif- 
ference with the SIS model is that once an individual gets 
infected it is removed and, therefore, it can not catch the 
infection again. If we start from a single infected node the 
size of the outbreak of the disease is given by the num- 
ber of nodes that can be reached assuming that each link 
is occupied with a probability A. Thus, the SIR model 
is equivalent to a bond percolation problem with bond 
occupation probability A |Q. Moreover, the size of the 
outbreak is just the size of the giant component. 

The study of the SIR model in random SF networks 
confirmed the epidemiological picture obtained for the SIS 
model in complex networks with power-law connectivity 
distributions [p5|,3(|. For instance, it has been shown that 
the effective epidemic threshold is inversely proportional 
to the connectivity fluctuations (fc 2 ) and hence it is van- 
ishing in the thermodynamic limit N — + oo for 7 < 3. 
The high heterogeneity of SF networks also causes that 
the relative incidence of an outbreak strongl y d epends on 
the connectivity of the first infected nodes |36|]. This de- 
pendency on the initial seed should not be confused with 
the dependency on the initial conditions found previously 
since for the SIR model the differences in the relative inci- 
dence of an epidemic outbreak is an intrinsic effect of the 
large heterogeneity in the connectivity distribution of SF 
networks and is not related to the presence of any topo- 
logical constrain. 

We have made numerical simulations of the bond per- 
colation problem on top of the structured networks gener- 
ated using the two variants of the deactivation model, A 
and B. Once the graphs are generated each link is removed 
with a probability 1 — A. Then, the size of the giant com- 
ponent of the resulting graph is computed. All the results 
reported below were obtained taking an average over 10 
graph realizations of the deactivation model and 10 real- 
izations of the link removal procedure. In Fig. || we plot 
the size of the giant component 5 as a function of A for 
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Fig. 3. Size of the giant component as a function of the bond 
occupation probability A for the network generated with m — 
3. The top figure (open symbols) corresponds to model A and 
the bottom figure (filled symbols) to model B. 



different graph sizes N for models A and B, with m = 3. 
Notice that in both cases the qualitative picture is the 
same. 

Between the limiting cases 5(0) = and 5(1) = 1, 
there is an intermediate range of A where 5 goes from a 
value close to to a value close to 1. The width of this 
interval decreases appreciable with increasing N. More- 
over, the point at which the transition takes place system- 
atically shifts to larger values of A approaching 1. Since 
the deactivation model has essentially a one dimensional 
topology we expect that in the large N limit 5 = for 
any A < 1. This hypothesis can not be confirmed by di- 
rect observation of Fig. || but it can be checked using a 
finite size scaling analysis. When we plot 5 as a function 
1 — A in a log-linear scale we observe equidistant shifts of 
the curves after increasing TV from 10 4 to 10 5 and then 
to 10 6 . This scenario corresponds to a scaling of the form 
5 = f[(l — X)N a ], where f(x) is a scaling function that is 
independent of N and a is a scaling exponent. 
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1-X 

Fig. 4. Size of the giant component as a function of the 
rescaled variable (1 — \)N a , with a = 0.15 and a = 0.14 
for models A and B, respectively. The symbols are in corre- 
spondence with those in Fig. ^| All the curves of model A were 
shifted by the same amount to the right to make evident the 
similarity between the curves of both models. 

In Fig. |J we plot S versus (1 — X)N a in a log- linear 
scale using the value of a that gives the best data col- 
lapse. This was achieved for a = 0.15(1) and a — 0.14(1) 
for models A and B, respectively. Within the numerical er- 
ror, the exponent a is the same for both models. Moreover 
the scaling function f{x) is also unique, up to a constant 
factor. These results point out that with regard to the 
bond percolation problem, or the SIR model, the differ- 
ences in the degree distributions of model A and B are 
qualitatively irrelevant. 



5 Discussion and conclusions 

We have study the dynamics of infectious diseases in struc- 
tured scale-free networks generated using the deactiva- 
tion model introduced by Klemm and Egufluz. By means 
of large scale numerical simulations, we have shown that 
the existence or not of an epidemic threshold for the SIS 
model depends on the initial density of infected individu- 
als when the connectivity fluctuations of the network are 
unbounded. This effect is completely new and reflects the 
peculiar topological nature of the networks generated with 
the deactivation model algorithm. Guided by the analyt- 
ical solution to the SIS model in an ensemble of discon- 
nected stars, we may summarize the SIS dynamics as fol- 
lows. 

The linear topological nature of the networks provoke 
the dynamics of the epidemic spreading to be almost purely 
diffusive, like in a one-dimensional lattice. However, for 
model B the small value of 7, that leads to the divergence 
of the connectivity fluctuations, allows the existence of 
stars with very large connectivities such that for any value 
of A the infection can always get trapped in these stars. 
Hence, starting from homogeneous initial conditions, the 



probability of hitting a star in which the epidemic survives 
is not zero. This is not certainly the case if one starts at a 
single infected node. In this case, there is no time for the 
infection to spread through the network. It would be of 
further interest to study this diffusion process in more de- 
tails, for example, by looking at the distribution of waiting 
times for each node in the cicle infected — > susceptible — > 
infected. On the other hand, when the connectivity fluc- 
tuations are finite (model A), the probability of finding 
a star with degree k > fc c (A) is zero below A c and thus 
the epidemic threshold is recovered for all possible initial 
conditions. 

Finally, the study of the SIR model confirms that the 
dominant factor in this case is the chain structure of the 
network with star-like nodes connected locally. The thresh- 
old for this epidemiological model coincides, as it should, 
with the critical point of the corresponding one-dimensional 
bond percolation problem that gives an epidemic thresh- 
old A c = 1. 

In summary, we have provided evidences that in this 
kind of networks the dominant factor determining the be- 
havior of both the SIS and the SIR models is the lin- 
ear topology of the network. The existence of high degree 
nodes introduces new effects such as the dependency on 
the initial conditions for the deactivation model B with 
7 < 3. As a final remark, we should say that the results 
for the SIS and the SIR models here reported cannot be 
directly extended to all correlated networks. The peculiar 
topological features of the networks generated with the 
deactivation model make them unique within the class 
of correlated networks. Some recent works have begun to 
[|5|[2|,[37j address the SIS and the SIR models in corre- 
lated networks establishing the general conditions for the 
existence or not of an epidemic threshold. The extensive 
numerical exploration of these models in such networks is 
still to do. 

We thank A. Vespignani and R. Pastor-Satorras for helpful 
comments and discussions. This work has been partially sup- 
ported by the European Commission - Fet Open project COSIN 
IST-2001-33555. 
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